Tunneling without tunneling: wavefunction reduction in a mesoscopic qubit 
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The transformation cycle and associated inequality are suggested for the basic demonstration of 
the wavefunction reduction in a mesoscopic qubit in measurements with quantum-limited detectors. 
Violation of the inequality would show directly that the qubit state changes in a way dictated by the 
probabilistic nature of the wavefunction and inconsistent with the dynamics of the Schrodinger equa- 
tion: the qubit tunnels through an infinitely large barrier. Estimates show that the transformation 
cycle is within the reach of current experiments with superconducting qubits. 
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Can a quantum particle tunnel through a barrier which 
has vanishing transparency? Immediate answer to this 
question is "no" as follows from the elementary proper- 
ties of the Schrodinger equation. More careful consider- 
ation should of course remind that evolution according 
to the Schrodinger equation is not the only way for a 
state of a quantum particle to change in time. Prob- 
abilistic nature of the wavefunction implies that it can 
also evolve due to the "wavefunction reduction" : ran- 
dom process of realization of one specific outcome of a 
measurement. In this process, the wavefunction of the 
system changes coherentlyfor any given outcome of the 
measurement - see, e.g, [H, 0, if measured with the 
mor ally best" or, in a more modern and descriptive 
language, quantum-limited detector. Such an evolution 
of the measured system can contradict the Schrodinger 
equation despite the fact that the dynamics of the mea- 
surement process as a whole is governed by this equation. 
It can be described formally as generic quantum oper- 
ation within the approach based on positive operator- 
valued measures (POVM) All "counter-intuitive" 
quantum-mechanical phenomena arise from the wave- 
function reduction. The best known example is given 
by the EPR correlations which violate the principle of 
"no action at-a-distance" as quantified by the Bell's in- 
equalities. They appear for specific random outcomes of 
the local spin measurements. On average, there is no 
action-at-a-distance, since the correlations do not vio- 
late the relativistic causality. In mesoscopic solid-state 
qubits, proposed and/or observed manifestations of the 
wavefunction reduction include violations of the tempo- 
ral Bell inequalities [1, @, 0, 0, measurements of the 



^10], stochastic reversal 
The purpose of this 



"weak values" of the operators 
of the wavefunction reduction 
work is to suggest a sequence of quantum transforma- 
tions and the corresponding inequality which illustrates 
directly the contradiction between the wavefunction re- 
duction and the Schrodinger dynamics by violating the 
very basic intuition that a particle can not tunnel through 
an infinitely large barrier. For mesoscopic structures with 
their small geometric dimensions, violation of this intu- 
ition would provide, arguably, more dramatic illustration 



of the wavefunction reduction than the non-locality of 
conventional Bell's inequalities. 

The system we consider is a qubit with the two basis 
states j = 0, 1, distinguished by the average values of 
some quantity x, for instance, electric charge or magnetic 
flux in the case of superconducting qubits. 

As the simplest example, x can be viewed as a position of 
an individual particle (electron in coupled quantum dots 
IgL irj , Cooper pair on a superconducting island 13, 13] , 
^HE quasiparticle in a system of two quantum antidots 
[13], or ultracold atom in a BEC junction 21]) 



which can be localized on the opposite sides of a tunnel 
barrier separating the states \j) and creating tunnel am- 
plitude A > between them (Fig. 1). The coordinate 
X is measured by a detector which converts the informa- 
tion about X into the classical output q. The detector is 
characterized by the probabilities Wj{q) of producing the 
output q, when the qubit is in the state For instance, 
in the example of the quantum-point-contact (QPC) de- 



tector (see, e.g., [22| and references therein), the qubit 



controls the scattering characteristics of electrons in the 
contact, and therefore the current / flowing through it 
(Fig. 1). In this case, the output q is the total charge 
transferred through the contact during the time r of the 
measurement. 
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FIG. 1: Illustration of the qubit measurement process, (a) 
Typical probability distributions Wj{q) of a detector output 
q for the qubit in the state \j), j = 0, 1. (b) Schematics of 
the particular measurement setup based on the QPC detector 
measuring mesoscopic "charge" qubit. In this case, the qubit 
states \j) differ by position of some elementary charge which 
controls the QPC transmission properties, and therefore, the 
current / in the contact driven by the applied voltage V. 
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The goal of the transformation cycle developed in this 
work is to demonstrate that the evolution of the wave- 
function in the measurement process is "real" to the same 
extent as the dynamics of wavefunction governed by the 
Schrodinger equation, i.e., it describes evolution of the 
physical quantities and not only information about them. 
This is accomplished by combining the two types of evo- 
lution in one cycle arranged so that they completely com- 
pensate each other, leaving the initial qubit state un- 
changed. The first part of the cycle is the wavefunc- 
tion reduction in a weak quantum-limited measurement 
which mimics the tunneling between the qubit states \j) 
through infinitely large barrier with A = 0. The sec- 
ond part is the regular tunneling with A 7^ 0. The fact 
that no charge or flux is transferred through the tun- 
neling barrier in the whole cycle means that the wave- 
function reduction induces tunneling even without the 
corresponding tunneling amplitude. 

The starting point of the cycle is a qubit with van- 
ishing average bias e = between the states The 
Hamiltonian of such a qubit is 



H = -{va, + Aa^)/2, 



(1) 



where ax,z are the Pauli matrices, and the bias v rep- 
resents the low-frequency noise characteristic for meso- 
scopic solid-state qubits 23, 24, [H, [2^. Due to as- 
sumed weak, but unavoidable, relaxation, the qubit at 
low-temperature T <C A is in the instantaneous ground 
state of ([T]) with the density matrix pi in the basis 

/ Co cqcA 
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where the probability amplitudes Cj are co,i = [(1 ± 
?;/f])/2]i/2, and Vl ^ {v^ + /S^f^. For A > ?;o, where 
vo is the r.m.s. magnitude of noise w, the state ([2]) is 
the ideal version of initial state for our transformation 
cycle, the eigenstate = 1 of the ax operator with 
co,i = l/v^. However, optimization of the cycle as a 
whole discussed below can require keeping the ratio A/wq 
finite, which reduces the noise-induced dephasing. 

The cycle begins by rapidly (on the scale h/vo) raising 
the tunnel barrier so that the tunnel amplitude vanishes, 
A — > 0, while the qubit wavefunction remains distributed 
between the states The next step is a weak measure- 
ment of the Gz operator performed on the qubit state 
(|2|) with a quantum-limited detector. Qualitatively, this 
means that the fiuctuations underlying the probability 
distributions Wj (q) of the detector output are themselves 
quantum, so that the evolution of the detector -|-qubit sys- 
tem leading to any given value q of the output is quantum 
coherent. An example of the quantum-limited detector is 
the QPC [22,l27| mentioned above (Fig. lb), in which the 
distribution of the transferred charge q is created by the 
quantum-coherent scattering of electrons at the contact. 
The information about the qubit state gained by the mea- 
surement depends on the observed output q (Fig. la), and 



implies that the qubit wavefunction ainplitudes Cj evolve 
into the g-dependent values Cj (q) (27l . [28| : 



C]\/wj{q) 



\/'i«o('?)|cop + wi{q)\ci\ 



(3) 



If the detector is not strictly quantum-limited, it in- 
duces partial dephasing of the qubit state even for a given 
specific detector output q. Physically, this dephasing is 
caused by the loss of information in the measurement 
process and can be described in general by the suppres- 
sion factor e^^ of the non-diagonal elements of the qubit 
density matrix p in the measurement basis. The qubit 
density matrix after weak cr^ measurement the is 



{II) Pn 



cUq), co(g)ci(g)e-''+»^W 
Ca{q)ci{q)e-^-^^<^-\ c^^{q\ 



(4) 



The factor e"^^'^^ here represents the noise-induced phase 
(^(t) = /g dtv{t)/h during the measurement time t. 

The aim of the next step of the cycle is to reverse the 
changes in the qubit amplitudes cj due to redistribution 
through the tunnel barrier with vanishing transparency 
in the reduction process ([3]) by regular tunneling. This is 
done by creating a non- vanishing tunneling amplitude for 
some appropriate period of time, i.e., realizing a fraction 
of the regular coherent oscillations in which the charge or 
flux goes back and forth between the qubit basis states. 
In the situation with no disturbances (vanishing noise 
and quantum-limited detector), this can be done pre- 
cisely, returning the qubit from the state 



dt 



\i^m) = WM<l)\0) + .ywi{q)\l)]/^Wo{q)+Wi{q). 

obtained as a result of the measurement, to the initial 
state = [|0)-|-|l)]/-\/2 before the measurement. Writ- 
ing the amplitudes of \ipm) as cos(0/2) and sin(0/2), one 
can see that the required transformation is the rotation 
about the y axis. For mesoscopic qubits, such a rota- 
tion corresponds to complex tunnel amplitude A' with 
argA' = 7r/2, and the rotation angle is 

Typically, the qubit structure allows only for the real 
tunnel amplitude A - see the Hamiltonian JT]), which re- 
alizes the X-axis rotations Rx = exp{i(Tx J A{t)dt/2h} of 
the qubit. In this case, the y-axis rotation ([5]), Ry — 
exp{—iay J \A'{t)\dt/2h}, can be simulated directly by 
the cc-rotation of the same magnitude, if it is preceded 
and followed by the z-axis rotations: Ry = R~^RxRz- 
The z-rotations Rf^ — exp{±i(T2 7r/4} are created by the 
pulses of the qubit bias: / e{t)dt/h — ±7r/2. Such a 
three-step sequence is simplified by interchanging the or- 
der and magnitude of rotations: first the x-, then one 
z-rotation: 



{III) 
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Under the experimentally-realistic assumption that the 
measurement takes much longer time than control pulses, 
the noise-induced distortions of the qubit states is dom- 
inated by the phase accumulation in Eq. ([¥]). The qubit 
density matrix after the pulses ^ then is 



(7) 



As the last step, one needs to check whether the qubit 
is brought back to the initial state Ux = 1. This is done 
directly by strong measurement of Ux, which gives, upon 
sufficiently large number of cycles, the error probability 
p of the qubit reaching the wrong state ax = —1: 



{IV) p = Tr[Pp\, P=\[\ / 



(8) 



If p = 0, the cycle (/) — {IV) returns the qubit with cer- 
tainty to the state tXa; = 1. This means that for any finite 
measurement strength, when the distributions wo{q) and 
wi{q) are not identical, the transfer of the qubit ampli- 
tudes due to the wavefunction reduction with A = is 
precisely compensated for by a fraction of a period of 
coherent qubit oscillations with A 7^ 0. Since the oscilla- 
tions actually transfer the charge or flux between the two 
qubit basis states, the fact that the cycle is closed shows 
that the qubit evolution in the wavefunction reduction 
can not be interpreted only as the changes in our knowl- 
edge of the qubit state, but rather involve actual transfer 
of charge or flux without the tunneling amplitude. 

In the presence of the finite noise and detector non- 
ideality, the error probability p is non- vanishing. Still, 
the probability 1 — p of the closed cycle can be larger 
than the value explainable by the classical description 
of the measurement process, demonstrating that the ac- 
tual qubit evolution in measurement is governed by the 
wavefunction reduction ^ . The alternative classical de- 
scription would be based on the assumption that the pro- 
cess of switching off the tunneling amplitude A does not 
leave the qubit in the state ^ but localizes it in one of 
the basis states on one or the other side of the tunnel 
barrier with some undetermined probability r, so that 
the qubit density matrix is p^"''^ = r|0)(0| + {l-r)\l){l\. 
The qubit state is then "objectively" well-defined and 
coincides with one of the basis states. It is, however, un- 
known, and the measurement provides information about 
this unknown state changing the probability r: 

r r{q) = rwo{q)/irwo{q) + (1 - r)wi{q)) . 

Applying the same transformations as to the state ([2]) 
to p^^^^q) = r{q)\0){0\ + (1 - r(g))|l)(l|, one finds the 
error probability p*-"^'-* {q) given the measurement outcome 
q. Averaging over the probability a{q) = rwo{q) -I- (1 — 
r)w2{q) of different g's, one obtains the total classical 
error probability 



P 



id) 



dq 



wa{q)wi{q) 
wo{q) + wi{q) 



(9) 



Equation ^ shows that the probability p*^"^'^ is inde- 
pendent of the assumed initial probability r, i.e. of 
noise v, and of the degree of detector non-ideality. The 
probability ([9|) characterizes the measurement strength, 
with p^"^'-* ~~> for strong projective measurements, when 
WQ{q)wi{q) = 0, i.e. the measurement provides definite 
information about the qubit state. For weak measure- 
ments, the two distributions wo.i{q) nearly coincide and 
pici) ^ 1/2. 

To see the incompatibility of the classical description 
with the actual qubit evolution, one needs to compare 

P' 



id) 



([9]) with the probability p found from Eqs. ([7]) and 

after averaging over the detector output: 



p = p(=') [1 - e~" F] 



(10) 



The factor F here describes effect of the noise v. 
F ' ^ 



( (.^ + A^)V2 cosy^ (^(Wr.)) . (11) 



Equation (I10|) shows that if the noise and the detector 
non-ideality are not very large: F —^ 1 and 77 — > 0, the 
errors of the quantum cycle are suppressed regardless of 
the measurement strength characterized by p^'^'' , and it 
should be possible to observe that 



p<p 



cl) 



(12) 



Demonstration of this inequality would mean that a 
larger number of the transformation cycles are closed 
than can be explained classically, without invoking the 
transfer of the amplitude Cj across the infinite barrier. 
Note that Eq. (fTO| is obtained for essentially arbitrary 
noise model, and is therefore valid in the general case of 
the quantum noise that corresponds to entanglement of 
the qubit ([1]) with an arbitrary quantum system averaged 
out in Eq. (fTTj) . 

To make the limitations on the noise more quantitative, 
we adopt the usual model of the low-frequency noise (see. 



e.g., [2J, l26[) as a classical Gaussian random variable 



static on the time scale of the cycle, which gives 



F : 



1 



2ttvq 



dve 



Acos{vT/h) 
(?;2 + A2)i/2 



(13) 



Equation p3)) can be evaluated analytically in several 
limits. For A ^ vo, the noise effect on the initial state 
disappears, and Eq. ()13|) reduces for arbitrary vot/H to 
pure dephasing 



(14) 



In the opposite limit of negligible dephasing, vqt/H ^ 1, 
or if dephasing is eliminated by a spin echo technique. 



F^ ^2a/ne''Koia) , a = A^/{4:vl), 



(15) 



where Kq is the modified Bessel function. This equation 
describes the transition from F ~ 1 for A ^ dq, to F ~ 



4 



(A/wo) ln(wo/A) for A < uq- For A < wq and t ^ 0, 
one can neglect the exponential factor in ([13]) to get: 



F^{2/7r)'^HA/vo)Ko{AT/h). 



(16) 



The strength of the qubit coupling both to the detec- 
tor and the noise depends typically on the same qubit 
parameters. This means that optimization for maxi- 
mum F should be done for fixed values of the vqt prod- 
uct which represents a more fundamental quality of the 
qubit-detector system. Figure [2] shows F as a function of 
the A/vq ratio for several values of vqt/H calculated nu- 
merically from Eq. ([TH)) . In the large- A limit, the curves 
saturate at F given by Eq. (|14p. which is exponentially 
small for strong noise. For vqt/K > 2, the A-dependence 
of F is non-monotonic, with the maximum at small A/vq 
representing suppression of dephasing by small A. Al- 
though the reduction of A increases the noise effect on the 
initial state, it effectively suppresses dephasing: the in- 
stances of stronger noise which would result in a stronger 
dephasing become irrelevant, since the qubit state is al- 
ready localized for them in one of the states \j) and does 
not contribute to the wavefunction reduction. On the 
other hand, the probability of obtaining realizations with 
small noise scales as A/uq, and F is also suppressed for 
A ^ 0. This gives the peak of F at small A/vq seen 
in Fig. [2] and for vqt/H 3> 1, described analytically by 
Eq. p6p . This peak may be useful in experimental real- 
ization of the cycle discussed in this work. 




FIG. 2: The noise suppression factor (|13|l as a function of the 
A/vo ratio for several values of the noise strength or qual- 
ity parameter vor/h of the qubit-detector system. For the 
discussion, see main text. 

In general, realization of this cycle and observation of 
the inequality is not very far from the possibilities 
of the current experiments with superconducting qubits. 
Indeed, Figure [2] shows that a noticeable difference be- 
tween the classical and quantum error probabilities can 
be obtained even for the cycle times r on the order of 
~ 3 dephasing times (|14p H/vq. Experimental dephas- 
ing times for the qubit states that differ by the average 
flux or charge values are about 5 ns - see, e.g., 1^ 29^ . 



giving roughly 15 ns interval for performing three opera- 
tions of the cycle when the qubit is subject to dephasing: 
measurement and two compensating pulses Simple 
qubit control pulses are regularly performed on the few- 
nanosecond time scale, and can be fit into this interval. 
The new requirement presented by the cycle described 
above is the need to perform a variant of the feed-back 
control, when the applied pulses depend on the result of 
the previous qubit measurement. Necessary rapid non- 
destructive read-out could be developed, e.g., using the 
superconductor qubit control circuits [28| . 

In summary, we have proposed a transformation cy- 
cle aimed at the very basic demonstration of the two 
main features of the wavefunction reduction in quantum 
measurements using mesoscopic solid-state qubits. The 
wavefunction can evolve in a way that explicitly contra- 
dicts the dynamics of the Schodinger equation: a parti- 
cle can be transferred through an infinitely large barrier. 
On the other hand, to the same extent as the Schodinger 
equation, the reduction affects not only the probability 
distributions of a dynamic variable (i.e., electric charge or 
magnetic flux), but the variable itself. Although the par- 
ticle transfer through the infinite barrier in this process 
is not identical to the regular Schodinger-equation tun- 
neling, e.g., there is no a corresponding current operator, 
the cycle developed in this work demonstrates that such 
a transfer can still take place and should be observable 
in current experiments with superconducting qubits. 

D.V.A. would hke to thank L.Y. Gorelik, A.N. Ko- 
rotkov, and V.K. Semenov for useful discussions. This 
work was supported in part by the NSA under ARO con- 
tract W911NF-06-1-217. 



[1] 
[2] 

[3] 



[4] 
[5] 

[6] 

[7] 

[8] 

[9] 

[10] 

[11] 
[12] 



G. Liiders, Ann. Phys. (Leipzig) 8, 322 (1951). 

M.L. Goldberger and K.M. Watson, Phys. Rev. 134, 

B919 (1964). 

J.S. Bell and M. Nauenberg, in: Preludes m Theoreti- 
cal Physics, (North HoUad, Amsterdam, 1966), p. 279; 
reprinted in: J.S. Bell, Speakable and Unspeakable m 
Quantum Mechanics (Cambridge Univ. Press, 1987). 
A. Peres, Quantum Theory, (Kluwer, 1993). 
A.J. Leggett and A. Garg, Phys. Rev. Lett. 54, 857 
(1985). 

A.N. Korotkov and D.V. Averin, Phys. Rev. B 64, 165310 
(2001). 

R. Ruskov, A.N. Korotkov, and A. Mizel, Phys. Rev. 
Lett. 96, 200404 (2006). 

A.N. Jordan, A.N. Korotkov, and M. Buttiker, Phys. 
Rev. Lett. 97, 026805 (2006). 

N.S. WiUiams and A.N. Jordan, Phys. Rev. Lett. 100, 
026804 (2008). 

A. Romito, Y. Gefen, and Y.M. Blanter, Phys. Rev. Lett. 
100, 056801 (2008). 

N. Katz et al, Phys. Rev. Lett. 101, 200401 (2008). 
Y. Nakamura, Yu.A. Pashkin, and J.S. Tsai, Nature 398, 
786 (1999). 



5 



[13] Yu.A. Pashkin et al, Nature 421, 823 (2003). 
[14] J.R. Friedman et al, Nature 406, 43 (2000). 
[15] J.H. Plantonborg et al, Nature 447, 836 (2007). 
[16] T. Hayashi et al, Phys. Rev. Lett. 91, 226804 (2003). 
[17] A.K. Huttel et al, Phys. Rev. B 72, 081310(R) (2005). 
[18] I.J. Maasilta and V.J. Goldman, Phys. Rev. Lett. 84, 
1776 (2000). 

[19] D.V. Avcrin and V.J. Goldman, Solid State Commun. 

121, 25 (2002). 
[20] S. Foiling et al. Nature 448, 1029 (2007). 
[21] D.V. Averin et al, Phys. Rev. A 78, 031601(R) (2008). 
[22] D.V. Averin and E.V. Sukhorukov, Phys. Rev. Lett. 95, 



126803 (2005). 
[23] G. Falci et al, Phys. Rev. Lett. 94, 167002 (2005). 
[24] G. Ithier et al, Phys. Rev. B 72, 134519 (2005). 
[25] R.H. Koch, D.P. DiVincenzo, and J. Clarke, Phys. Rev. 

Lett. 98, 267003 (2007). 
[26] R. Harris et al, Phys. Rev. Lett. 101, 117003 (2008). 
[27] A.N. Korotkov, Phys. Rev. B 60, 5737 (1999). 
[28] D.V. Averin, K. Rabcnstein, and V.K. Semenov, Phys. 

Rev. B 73, 094504 (2006). 
[29] T. Duty et al, Phys. Rev. B 69, 140503(R) (2004). 



